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Abstract 

We describe a simple mechanical system that involves Spontaneous Symmetry Breaking. The system consists 
of two beads constrained to slide along a hoop and attached each other through a spring. When the hoop rotates 
about a fixed axis, the spring-beads system will change its equilibrium position as a function of the angular 
velocity. The system shows two different regions of symmetry separated by a critical point analogous to a second 
order transition. The competitive balance between the rotational diynamics and the interaction of the spring causes 
an Spontaneous Symmetry Breaking just as the balance between temperature and the spin interaction causes a 
transition in a ferromagnetic system. In addition, the gravitational potential act as an external force that causes 
explicit symmetry breaking and a feature of first-order transition. Near the transition point, the system exhibits a 
universal critical behavior where the changes of the parameter of order is described by the critical exponent [3—1/2 
and the susceptibility by 7 = 1. We also found a chaotic behavior near the critical point. Through a demostrative 
device we perform some qualitative observations that describe important features of the system. 



I Introduction 

Although gauge theories and Spontaneous Symmetry 
Breaking (SSB) are rather abstract concepts which are 
associated with symmetries operating on internal degrees 
of freedom, their significance and nature emerges in a quite 
straightforward way through systems that exhibit geomet- 
rical properties associated with space-time symetries, and 
that shares the same basic features as models with gauge 
structure. In the literature some works about analogies 
with classical systems are found 0, where different me- 
chanical models are used to illustrate similar characteris- 
tics found in particle physics and solid state physics. In 
this work with a simple mathematical treatment, we de- 
scribe another mechanical system that exhibits SSB and 
dynamical restitution of a discrete symmetry in a similar 
fashion as a second order phase transition. In contrast to 
the systems reported by the literature, our model is in a 
closer relation with a thermodynamic system in the sense 
that the former use the gravity as an essential external 
parameter to generate symmetry restitution, while ours 
may generate spontaneous breaking and restitution with 
only internal interactions, where the gravity act as an ex- 
ternal field that causes explicit symmetry breaking, such 
as external interactions cause explicit breaking in thermal 
physics. The model shows not only interesting features 
that resembles to the SSB in gauge theories, but also in- 



volves other related concepts. For example, it can be used 
in a basic level of classical mechanics as an application of 
rotational dynamics and to visualize the stability of sys- 
tems from energy diagrams. In a higher level, it can be 
introduced as an example in the Lagrangian formulation 
of analytical mechanics. In a thermodynamic framework, 
this model illustrates the behavior of a first- and second- 
order trasition according to the Tisza's definition. Futher- 
more, because of the nature of the SSB, this system ex- 
hibits a critical behavior near the transition point which 
belongs to the same universality class as a large number 
of thermodynamic systems. Thus, this problem may pro- 
vide instructive insights that concern with more advanced 
topics into the classical study of the critical phenomena 
in phase transitions. In addition, the system presents in- 
teresting properties associated with nonlinear solutions, 
which may be source of chaotic behavior and generates a 
simple application into the theory of classical chaos. 

II Stability 

We consider two identical beads of mass m that can slide 
without friction on a horizontal hoop of radius R. There 
appears a potential of interaction through a spring at- 
tached to both beads. The spring constant is k and its 
equilibrium length is less that the diameter of the hoop, 
i.e. 2ro < 2R. The system is sketched in Fig. where the 
position zq of both beads is the same measured from the 
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Figure 1: BHS model. Two beads attached by a spring slide with- 
out friction along a horizontal hoop with rotation. The beads have 
the same position zq from the center. 



center of the hoop. There are two equivalent equilibrium 
positions determined by the transformation zq — > — zo, 
which illustrates the presence of a discrete symmetry. The 
hoop rotates about the z axis at angular velocity which 
causes the equilibrium position zo to decrease and the 
spring to stretch from its natural length 2ro to a longer 
length 2r. The lagrangian contains the kinetic energy of 
each bead and the elastic potential of the spring, while 
the gravitational potential is not considered because the 
center of mass of the system lies along the horizontal axis 
z. It is more suitable to write the lagrangian in cylindrical 
coordinates, which displays in a simpler form the symme- 
tries and constraints of the system. The lagrangian is then 
written as 



m(r + r <j> + z ) — 2k(r — ro) 



(i) 



The change in time of the azimuthal angle <j) is identi- 
fied as the angular velocity of the hoop. We can eliminate 
the coordinate r through the constraint r = ±\J R 2 — z 2 
which indicates that the motion of the beads is restricted 
along the hoop. Hence, the lagrangian can be expressed 
in terms of the generalized coordinate z as 



1 .2 
-\iz 



Vef(z), 



(2) 



where we have defined an effective mass [i — 
and a one-dimensional effective potential given by 
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Figure 2: Effective potential as a function of the z coordinate for 
angular velocity (A) lj > uj c , (B) uj = uj c and (C) uj < uj c - 
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To investigate the stability of the solutions above, we 
must examine the second derivative of the effective poten- 
tial, which gives 



2 V e 



ef 



dz 2 



■ 4kr 



R 2 



(R 2 - z 2 ) 



2^3/2 



(6) 



This expresion evaluated in the 3 equilibrium points, 
displays different possibilities according to the sign of the 
second derivative. The results are resumed in table |IJ 
where R 2 > £ 2 was d emanded in order to assure real so- 
lutions for z = ±^R 2 -^ 2 . In Fig. we show three 
different plots of the effective potential as a function of 
the coordinate z, where we distinguish a critical angular 
velocity determined by 



2 2k(R-r ) 

UJ r = 



rnR 



(7) 



and two regions that depend on the value of the angular 
velocity of the hoop. For values with uj 2 > uj 2 (symmetri- 
cal phase), there is just one minimum at z = that corre- 
sponds to the case in which both beads hold their maxima 



V e f(z) = Vel ~ V ce nt = 2k ( ^ R 2 - Z 2 - T ) - 7TIUJ 2 (R 2 - Z*) , 

(3) 

with V e i the elastic potential of the spring and V cen t 
the centrifugal potential which contains the rotational dy- 
namics. It is evident that the lagrangian in (0) does not 
change under the coordinate transformation z — » — z. The 
introduction of the effective potential in Eq. © will help 
us to determine the stability behavior of the system. Equi- 
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Table I: Stability regions according to the angular velocity. There 
is a critical angular velocity that separates two different phases. 
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Figure 3: Analitycal behavior of the effective potential with the 
angular velocity, (a) and (b) show the continuity of V e f and —^f- • (c) 
displays an infinite discontinuity across the critical point uj c (in these 
plots is 23.26mci/-s), exhibiting a second-order phase transition. 



distance in a symmetrical position. When the angular ve- 
locity slowly decreases near to the value determined by 
uj 2 — oj 2 , the minimum of energy becomes flatter. Im- 
mediately below the critical value uo Cl the single minimum 
bifurcates into two degenerated minima z = 
with the same energy. However, the spring chooses one of 
the minima causing an spontaneously symmetry breaking, 
where the symmetry transformation z — » — z is not man- 
ifest, and the spring comes to a new equilibrium position 
in either a positive or a negative coordinate z, which is 
greater with decreasing oj. Thus, the absolute value of z 
measures the degree of breaking just like an order param- 
eter. 

We can describe this process like a second order phase 
transition if we investigate the behavior of the effective 
potential in the critical point with respect to the variation 
of the angular velocity. The Ehrenfest definition estab- 
lish that a phase transition with respect to a thermody- 
namic variable (like temperature) is of n th order if the 
energy and its (n — 1) st-order derivatives are continuous 
at the transition point, whereas the n th -order derivative 
suffers a finite discontinuous jump. In fact, our system 
presents continuous derivative at first-order but suffers an 
infinite discontinuity in the second-order derivative d ^j/ 
as is shown in Fig. 03 where the angular velocity oj works 
like the temperature in a thermodynamic system. The in- 
finite jumps are described by the Tisza's Theory j^j which 
describes most of the higher-order phase transition. 

We emphasize that this behavior is analogous but not 
equal to the SSB in quantum theory. The analogy arises 
when the minimum of energy splits into new degenerated 
minima. Nonetheless, classically the system will maintain 
its symmetrical unstable equilibrium unless there appears 
external fluctuations that are not being considered into 
the equation of motion (in this sense this is not a "true" 
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Figure 4: Effective potential as a function of the z coordinate for 
angular velocity (A) uj > uj c , (B) uj = uj c and (C) uj < uj c . The 
gravity causes the negative minimum to be lower. 
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Figure 5: Projections of phase diagram in the (a) zq — uj and 
(b) zo — Sina planes. The " iso- velocities" in (b) show a first-order 
transition for uj < uj c , where zq jumps discontinuously between A 
and B. 

SSB). On the other hand, the quantum theory introduces 
fluctuations associated with the uncertainty principle in- 
herent into the structure of the theory. Such fluctuations 
appear spontaneously to unbalance the system toward a 
not symmetrical state . 

If the horizontal hoop is turned into an angle a, the 
gravity induces an explicit symmetry breaking. In this 
case the effective potential in changes as 



V e f(z) = 2/c(\AR 2 



ro 



(R 2 — z 2 ) + 2mgzSma, 
(8) 



where the last term corresponds to the gravitational 
potential of both beads. The zero gravitational poten- 
tial is taken in the symmetrical point z = 0. It is clear 
that Eq. does not hold the symmetry z — > — z due 
to the gravitational term, which will favor negative po- 
sitions with lower energy. In this case, the symmetrical 
plots from Fig. [2] will suffer slight deformations, such as 
is illustrated in Fig. ^ where we can see that for uj < uj c 
the minima splits into two different values as local and 
global minima. This behavior is observed, for instance, 
in an atom or molecule that presents a degenerated spec- 
trum of energy due to the spherical symmetry. If there 
are external fields such as electric or magnetic fields, this 
spectrum splits into many distinguished levels of energies 
(Stark and Zeeman effect). We can obtain a phase diagram 
from the minima conditions. When the axis of rotation is 
in horizontal position, the minima solutions in Eq. ® as 
a function of uj describes a plot analogous to the 
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magnetization-temperature phase diagram in a ferromag- 
net, such as Fig. Efa) shows. If the z-axis holds different 
angles a with respect to the horizontal plane, the min- 
imum condition from the effective potential in Eq. ((SJ 
describes a diagram as a function of Sina similar to the 
isotherms in a magnetization-magnetic field phase dia- 
gram (see Fig. EJb)). When uo < uo Cl Fig. I^b) illustrates 
how the global minimum jumps discontinuously from — zq 
to zo when a passes through zero, like a first-order tran- 
sition. 

Ill Critical Exponents 

Above, we identified the equilibrium position zq as the 
parameter of order which measures the degree of breaking 
of the symmetry. The behavior of second-order transition 
suggests us a critical behavior near the transition point. 
The effective potential from Eq. Q can be written as 

Vef(z) 2k (R 2 + rl) - muo 2 R 2 + (muo 2 - 2k) z 2 



f/-'n/M I - ( ^J • (9) 



2k(R-r ) 



Since nearly below the critical point uo c — mR 
the minimum at zo = suffers small deviation z ~ ±&z, 

we can expand the root term in as \J 1 — (j=>) 2 = 

1 — | (^) 2 — \ (^) 4 — ' ' "> so that the effective poten- 
tial takes the form of a Landau-type expansion near the 
critical point 

Vef(zo) = do + CL 2 Z 2 + d±Z^ ~\ , 



where the coefficients up to fourth order are 
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Furthermore, near uo c the coefficient in (|llb|) becomes 
<i2 ~ 2mcj c (cj — o; c ) = a° (a; — a; c ) • (12) 

The Landau theory ^ is formulated in the framework 
of the mean-field approximation, which in our case takes 
the form of a minimum condition. Thus, the potential in 
(fTU|) up to fourth order must accomplish the condition 



dVef(zo) 

dz 



Zo [2(12 (uO — UO c ) + 4<24Zo] = 0- 



(13) 



For uo > uo c the only real solution is zq =0, while for 
uo < uo c two real minima solutions are obtained 



zo 



2a4 



1/2 



(14) 



Thus, the parameter of order spontaneously becomes 
nonzero and grows as ^uo c — uo for angular velocities nearly 
below uo c . The critical exponent (3 = - associated with 
the order parameter is identical to those obtained by the 
classical Landau theory for thermodynamic parameters 

Ej, such as the magnetic moment in a ferromagnetic 
substance, the difference in Zn-Cu occupation in a binary 
alloy, the molar volume liquid-gas in a fluid, among other 
second-order transitions. 

This system could be thought as a mechanical equiva- 
lent of a thermodynamic system, where the state is char- 
acterized by three "thermodynamic" variables: the equi- 
librium position zq (parameter of order), the angular ve- 
locity uo (temperature-like variable) and the gravitational 
force (pressure-like variable). A change of state can take 
place if we introduce a " heat reservoir" , which in our case 
is a power source (for example an electric motor) that 
provides the rotational energy to the hoop. If we con- 
sider the complete potential including the gravitational 
coupling W = 2mgSina, the condition in Eq. (p"3|) takes 
the form 

dVef d [ Z ^ = 2a§ (uo - uo c ) zo + 4a 4 ^o + W = 0. (15) 
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Table II: Comparison between mechanical and ferromagnetic systems. M is the magnetic moment, kB the Boltzmann's constant, 
jjio the Bohr magneton and q the coupling constant of a mean field (which is the magnetic equivalent to the spring constant k) . 
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Fi glire 6 1 Frecuency of oscillation as function of uj. At uj c , the 
frecuency becomes zero. Above uj c the frecuency f2> grows with uj 
without limit, while below, the frecuency f2< grows as uj decreases 

The mechanical susceptibility, which is understood as the 
linear response of the parameter of order (minimum po- 
sition zq) due to an infinitesimal conjugate field (W) is 
Xmech — ff^- This can be calculated if we differentiate 
the Eq. (|15|) with respect to W, obtaining 



Xmech 



-1 



2a\ (oj — ujc) + 12a4^Q 



for U! > U! c 

for uj < lo c 



(16) 



from which we deduce the critical exponent 7 = 1. Ta- 
ble m lists several corresponding features between the me- 
chanical and the ferromagnetic system in the framework 
of the two-dimensional Ising model under the Mean-Field 
approximation 6 . The results are listed in terms of the 
equilibrium positions at uj = (z (u;^o)) an d the magne- 
tization at T = (M( T ^ ))- 

IV Small Oscillations 

We can study small oscillations about the equilibrium po- 
sitions. For the symmetrical phase with uj > uj c , the effec- 
tive potential in (JHJ) can be expanded in a Taylor's series 
about the minimum at z = 0, from where we obtain an 
harmonic potential given by 



Vef(Sz) w cons. + 2kr ( - i 



(Sz) 



: cons. H — — [oz) . 



(17) 

where we can identify an effective spring constant k e f = 
4/cro — . We can even identify the effective mass de- 
fined by the Eq. @ at z = as \i = 2m. Then, the 
frequency of oscillation above the critical point is Vt > = 

ypf = yj 2 ^ where 1^1 > \ R \- The frecuenc y 

Q> approaches zero near uj c1 and the period of oscillation 
becomes infinite in the threshold of stability at z = (see 
Fig. EJ. This frecuency increases without limit with uj. 
For the breaking phase with uj < uj Ci we ca n perform a 
similar analysis about the minima at z = R 2 — £ 2 . In 

this case, the effective constant is fc e j = 4/cro — |^ 
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Fi gure 7: Graphical solutions for < cj c - The position as a func- 
tion of time is shown for (a) small perturbations at z > 0, greater 
perturbations just (b) below and (c) above the jump over z = 0, 
and (d) for enormous perturbations. Each solution (A,B,C and D 
respectively) is shown (e) in the phase space. 

and the effective mass is /i = 2m jr. Hence, the frequency 



is Q< 



— I 2kr 



^~ — -^2^ where |£| < \R\ . Here, as uj de- 
creases, the function decreases and reaches its mini- 
mum value at £ = ro when the hoop comes to a complete 
rest with uj = 0. As the angular velocity decreases, the fre- 
cuency Q< below the critical value becomes greater and 



reaches its maximum value ft < = 
uj = 0. 

V Chaotic Behavior 
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The complete description of the system is obtained as so- 
lution from the Lagrange's equations associated with the 
generalized coordinate z 

W 9L =0. (18) 
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With the lagrangian defined by Eq. (0) and the ef- 
fective potential from Eq. J3J), we obtain the equation of 
motion 
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Figure 8: Graphical solutions for uj > uj c . The position as a func- 
tion of time is shown in (a) for (A) small perturbations and (B) 
greater perturbations. Each solution (A and B respectively) is shown 
(b) in the phase space. 



which is a quite complex equation to solve. We note 
that this equation involves non-linear terms, which sug- 
gests us a chaotic behavior. From data given in Sec. IVIl 
we consider numerical solutions with different initial po- 
sitions around the equilibrium points (the entrance z(0) 
was always taken zero). Figs, in display the behavior of 
the position z as function of time and in the phase space 
(z Vs. z) for uo < uo c . We see clearly that for small per- 
turbations, the oscillations are restricted in regions with 
z > 0. In the phase space, this motion describes elliptic- 
like paths (we also see that these paths are not perfect 
ellipses due to anharmonic terms). If we increase the per- 
turbation, the system reaches the transition point, where 
the spring just oversteps the local maximum at z = 0, 
and the oscillations extend to both regions. For greater 
perturbations, the motion increases around all the hoop. 

On the other hand, we can see in Fig. [5| that for 
uo > uo Cl the motion is developed around the equilibrium 
point at z = 0, in agreement with the results of Sec. IIVI 
Figs, in [U| represent solutions for a value below but near 
the critical value, where we introduced small perturba- 
tions. For a slight perturbation (dzo ~ 5.92 x 10 _3 m), the 
system develops a motion in the breaking region. But if 
this perturbation is changed into Szo ~ 6.926 x 10 _3 m, the 
system presents random oscillations between both regions, 
where the graphs in the phase space describes chaotic 
paths. The pattern of this chaotic behavior is in extrem 
sensitive to very tiny changes of initial conditions (with 
differences of just ~ 1 x 10 _13 m). Exactly on the crit- 
ical value (Fig. ITU)) , the system presents very long and 
irregular oscillations. 

VI Experimental Description 

We made a demostrative BHS-apparatus using a hoop 
with radius R = 6 cm, a soft spring with natural length 
ro = 4.6 cm and constant k = 8.81 N/m, and two iron 
rings used as beads, each one with mass m = 7.6 g. The 
hoop was coupled to the axis of a DC motor with ad- 
justable velocity (which will work like a "heat reservoir"). 

Aligning the rotational axis in a horizontal position, 
the system describes the two regimes of symmetry. The 
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Figure 9: Graphical solutions for uj slightly below u) c . The position 
as a function of time is shown for (a) very small perturbations at 
z > and (c),(e) greater perturbations. Each solution ((b), (d) and 
(f) respectively) is shown in the phase space. These plots exhibit a 
chaotic behavior. 

experimental observation is illustrated in Fig. El accord- 
ing to the angular velocity. The measured critical velocity 
was uo c = 25.24 ± 2.1 rad/s (241 ± 20 RPM), which is 
in a quite good agreement with the results from table [IJ 
where a critical point at uo c = 23.26 rad/s (222.1 RPM) 
is expected. In addition to the motion at uo c1 Fig. ITTI 
shows two cases with uo < uo c and one case with uo > uo c . 
When the velocity was decreased below the critical value, 
we noted that the chosen direction {-\-z or — z) was highly 
sensitive to the way that the velocity was reduced and to 
the mechanical fluctuations near the transition. When the 
velocity was slowly reduced in a quasi-static process across 
the transition, the equilibrium position was maintained at 
z = (although we lubricate the hoop, the unstable equi- 
librium found some stability due to the small friction), but 
after some hesitating seconds the system finally fell into 
one of the stable minima (long relaxation time), break- 
ing the symmetry. However after several trials, we could 
see that the symmetry was hidden in the sense that we 
obtained about the same number of breakings into +z as 
into —z, i.e. both posibilities have the same probability of 
occurrence. Thus, we can see how the symmetry breaking 
is spontaneous, where the non-symmetrical states present 
two equally probable possibilities with the same energy. 
On the other hand, the device had a high sensitivity to 
the precision of the horizontal alignment. When the ro- 
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tational axis was slightly deviated with uo > uo c , the sys- 
tem always shows a breaking in one direction due to the 
gravitational potential, in agreement with the behavior 
described by Fig. 21 We estimated experimental explicit 
breaking within a range a < 0.1 rad (< ±5°). 
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Fi gure 10 1 Graphical solutions for uj = uj c . The position as a func- 
tion of time is shown for (a),(c) small perturbations. Each solution 
((b) and (d) respectively) is shown in the phase space. These plots 
exhibit a long oscillation with chaotic behavior. 



VII Summary and Conclusions 

In this work we studied a bead-hoop-spring (BHS) system 
under rotation as was described through the sections. The 
system presents different stability situations according to 
the angular velocity of the hoop, exhibiting a second or- 
der transition behavior, and where the position zq of the 
center of mass measured from the center of the hoop de- 
fines a parameter of order. The symmetry transformation 
z — > — z of the lagrangian is manifest if the angular ve- 
locity is above a critical value uo c1 where the equilibrium 
position at zq = is stable. Below uo Cl this equilibium 
becomes unstable, and appears two degenerated minima 
with zo 7^ 0. Thus, with only internal interactions (spring 
force and normal force between hoop and beads), the sys- 
tem exhibits an SSB below the critical value. If the gravi- 
tational potential is considered, the stability regions suffer 
a deformation, changing the equilibrium position at z = 
into another different from zero, and split ing the minima of 
energy at zq ^ into two non-equivalent minima. Under 
these circumstances, the system shows a preferential di- 
rection exhibiting explicit symmetry breaking and a first- 
order transition behavior. The experimental observations 
were in quite good agreement with the theoretical descrip- 
tion, where the critical velocity was measured within an 
8% of precision. 

The system belongs to the same universality class that 
many thermodynamic systems, exhibiting a Landau ex- 




Fi gure 11: Experimental device for frecuency of rotation (a) / = 
178RPM (a = 18.64md/s), (b) / = 210RPM (u = 21.99rad/s), 
(c) / = 246RPM (uj = 25.13rad/s « w c ) and (d) / = 366RPM 
(<jj = 38.32rad/ s > ex- 



pansion. From a mean-field-like approximation, we ob- 
tained the critical exponent (3 = \ and 7 = 1 associated 
with the parameter of order and the "mechanical suscep- 
tibility" respectively. 

The frecuency of small oscillations about the stable 
minima was also calculated. Above uo c , the frecuency Vt > 
increases without limit with uo. Near the critical value, 
the motion about z = describes long-period oscillations, 
which becomes into unbounded motion when uo passes 
through the critical point. Below uo c , the frecuency ^< 
about the new minima increases with uo decreases. 

The equation of motion exhibits some chaotic behav- 
ior near the critical value, where the system reacts in a 
random way for similar perturbation entrance. 

This analogy between the mechanical system and a 
thermodynamic system is highly suggestive, and it is very 
tempter to apply a deeper study of fluctuations below the 
critical point through the definition of a correlation func- 
tion and the implementation of a fluctuation-dissipation 
theorem [5] , which could provide some insights about 
the physical sense of these concepts. 
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